NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



1MSC-HRIC TR 0784480 

\S A-Cii- 168b 10) MANUFACTURING IN SPACI: N82-19236 

JIB DYNAMICS NU MERIC AL ANALYSIS Annual 
> .port, Aug. 1980 - Aug. 1981 (Lockheed 

li ssiles and Space Co.) 51 y HC A04/MF A01 Unclas 

C SC L 22A GJ/12 15118 


MANUFACTURING IN SPACE: 
FLUID DYNAMICS NUMERICAL 
ANALYSIS 


Augvsf 198 i 


Correct HASW-3281 (Annua! Report) 


Prepared for 

NASA HEADQUARTERS 
WASHINGTON, DC 20546 


. ■>» 

5. i. Robertson 
L A. Nicholson 
L W. Spradley 



Lockheed Missiles l Space Company, Inc. 
Huntsville Research A Engineering Center 
4800 Bradford Drive, Huntsville, Al 35807 


L.MSC-HREG TR 0 / B-H BO 


FOREWORD 

This document is an annual report describing the results of effort by 
personnel of Lockheed Missiles & Space Company, Inc., Huntsville Research 
& Engineering Center, for the National Aeronautics and Space Administration 
under Contract NASW-3281, u Manufacturing in Space: Fluid Dynamics Nu- 
merical Analysis," The contractual effort described in this document was 
performed during the year from August 1980 to August 1981 , The NASA 
Technical Director for this contract is Dr, Robert F. Dressier, Manager, 
Advanced Technology Program, NASA Headquarters, Washington, DC. 
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INTRODUCTION 


There is considerable current interest in the possibility of utilizing the 
near-zero gravity environment aboard orbiting space stations for various 
materials processing applications. There are a number of potential appli- 
cations where the near absence of convective stirring during the process 
would presumably enhance the quality of the product. The purpose of the 
research program described in this document is to perform numerical com- 
putations of thermally induced convection for various experiment configura- 
tions under microgravity conditions simulating the orbiting space station 
environment. The research program was initiated in August 1979. This 
document is an annual report summarizing the results of effort during the 
year from August 1980 to August 1981 . 


Three basic tasks were performed during this year 1 s effort: 

1. An analysis was performed of thermal convection in the Lal/Kroes 
experiment configuration. The analysis was based on a simplified 
geometry where three-dimensional flow conditions in the cubical 
shaped container were approximated by axisymmetric flow in a 
spherical shaped container of equal volume. 

2. An analysis was performed of thermal convection in a two- 
dimensional circular enclosure where the gravity vector changes 
in both magnitude and direction. 

3. An analysis was performed of the effect of container shape on the 
intensity ox natural convection. Results were obtained for two- 
dimensional enclosures of circular, half-circular and square 
shaped cross sections. 


The results of these tasks are described in the following Parts 1, 2 and 
3 of this document. 
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Part 1 

NUMERICAL SIMULATION OF NATURAL CON- 
VECTION IN A SPHERICAL CONTAINER DUE 
TO COOLING AT THE CENTER (IDEALIZATION 
OF THE LAL/KROES EXPERIMENT) 
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ABSTRACT 

Natural convection in a spherical container with cooling at the center 
was numerically simulated using the Lockheed-developed General Inter polants 
Method (GIM) numerical fluid dynamics computer program. The numerical 
analysis was simplified by assuming axisymmetric flow in the spherical con- 
tainer, with the symmetry axis being a sphere diagonal parallel to the gravity 
vector. This axisymmetric spherical geometry was intended as an idealiza- 
tion of the proposed Lal/Kroes crystal growing experiment to be performed 
on board Spacelab. Results were obtained for a range of Rayleigh numbers 
from 25 to 10,000. For a temperature difference of 10 C from the cooling 

sting at the center to the container surface, and a gravitational loading of 

-6 , -5 / 

10 g e# a computed maximum, fluid velocity of about 2.4 x 10 cm/sec was 

reached after about 250 sec. The computed velocities were found to be ap- 
proximately proportional to the Rayleigh number over the range of Rayleigh 
numbers investigated. 


m 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D784480 


CONTENTS 

Page 

FOREWORD 

ABSTRACT Hi 

NOMENCLATURE v 

INTRODUCTION 1 

NUMERICAL SIMULATION 2 

RESULTS 5 

CONCLUSIONS 14 

REFERENCE 15 


Lv 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D784480 


Symbol 


NOMENCLATURE 

Description 


g 

Se 

r 

R 

Ra 

T 

AT 

t 


t 

v 

V 

a 

P 

v 


gravity force in units of acceleration 

2 

gravitational acceleration on surface of earth, 980 cm/sec 
radial distance 
sphere radius 

3 

Rayleigh number = ^ ^ 

temperature 

temperature difference between cooled crystal growing 
surface at center of container and container walL 

time 

2 

dimensionless time = pt/R 
velocity 

dimensionless velocity j v 

gp AT R^ 

thermal diffusivity 

volumetric coefficient of thermal expansion 
kinematic viscosity 
stream function 
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INTRODUCTION 

One of the proposed experiments to be performed using the Fluids Expe*!' 
ment System (FES) facility aboard Spacelab Is the Lal/Kroes experiment (Ref. 1), 
which Is concerned with the controlled growth of a single crystal under near- 
zero gravity conditions. The crystal Is to be grown from a trlglycine sulfate 
(TGS) solution contained In a 10 x 10 x 10 cm cubical shaped container. The 
concentration Is to be nominally 45 g TGS/lOO cc water. The temperature 
of the solution will be within the range of 35 to 50 C. The crystal will be 
grown on a cooled 1 utn diameter disk positioned at the center of the container 
at the end of a 2 cm diameter insulated sting. The disk will be cooled up to 
10 C less than the surrounding container walls. 

The purpose of the experiment is to study crystal growth in the absence 
of significant gravity induced convective stirring of the solution. The near 
elimination of convective stirring is expected to result in the growth of high 
quality crystals. The purpose of this numerical study is to predict the in- 
tensity of convective stirring due to the small residual gravity forces re- 
maining under orbital flight conditions. The Lockheed developed General 
Interpoiant Method (GIM) code computer program (Ref. 2) was used in the 
numerical simulations. This computer code numerically integrates the basic 
fluid dynamics equations in conservation form. Computations were performed 
on the NASA- Langley Cyber 203, 
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NUMERICAL SIMULATION 

To permit reasonable economy in computer usage, the cubical shaped 
experiment configuration was modeled by a sphere of the same volume as 
the cubical container, with axisymmetric flow assumed. The axis of sym- 
metry is along a sphere diameter parallel to the gravity vector. The nodal 
point distribution was generated by the GIM code geometry module by speci- 
fying an array of 20 x 10 area elements and 21 x 11 nodal points over the 
half-sphere cross section enclosed by a semicircular arc and sphere axis 
of symmetry. The geometry module treated the semicircular region as a 
four-sided figure, with the bottom side the axis of symmetry and the other 
three sides concentric circiv ; ar arcs. The region was divided into an array 
of quadrilateral elements with curvilinear sides, and with the nodal points 
located at the corners of the individual elements. The geometry is shown 
in Fig. 1 with the computational grid network superimposed. The cooled 
crystal growing surface is represented by the outline of grid lines at the 
center of the sphere. The gravity vector is directed downward as shown 
in the figure, 

A steady state conduction temperature distribution is generated by the 
computer program as an initial condition. The nodal points outlining the 
cooled crystal growing surface are maintained at a temperature 10 C less 
than the surrounding container surface. 

The physical properties of the TGS solution were not well known at the 
time of this analysis. Based on conversations with Dr. R. L. Kroes, the NASA 
Principal Investigator, physical property values were assumed as listed in 
Table 1. 
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Fig. 1 


- Geometry and Computational Grid for Lal/Kroes Experiment 
Numerical Simulation 
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Table 1 

LIST OF ASSUMED PHYSICAL PROPERTY VALUES FOR TRIGLYCINE 
SULPHATE (TGS) SOLUTION (45 g TGS/lOO cc WATER) 


Viscosity, n 

1,78 centlpolse 

Thermal Conductivity, k 

0,00143 cal/cm-sec-C 

Density, p 

1.15 g/cm^ 

Specific Heat, C^ 

KBSESSSBSBQH^ 

Thermal Expansion Coefficient, 
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RESULTS 

Numerical simulations were performed for Rayleigh numbers, Ra, vary- 
ing from 25 up to 10,000. The computed spatial maximum velocities are plotted 
in JTig. 2 as a function of time for the various Rayleigh numbers. Note that by 
nondimens ionalizing the velocity and time as shown, the results collapse very 
closely about a single curve for Rayleigh numbers through 2500, The Ra r 
10,000 results are in accord with the other Rayleigh number results through 
the initial transient up to a dimensionless time of about 0.1, 

For the lower Rayleigh numbers, the results shown In Fig', 2 indicate the 
following approximate relationship for the spatial maximum velocity as a func- 
tion of time* 

v max * °- 00485 4>] 

< 0.00485 Ra [l - exp (-42 -^)J 

For the Lal/Kroea experiment configuration, the equal volume sphere radius, 

R, corresponding to the 10 x 10 x 10 cm cube is 6,2 cm, The time required to 
reach steady state is about 250 sec after imposition of an impulse gravitational 
load, The maximum spatial velocity at steady state is 2.4 x 10”® cm/sec for 
a 10"^ g gravitational load and a 10 C AT. 

The temperature, streamline and velocity contours at steady state are 
shown in Figs. 3 through 5 for various Rayleigh numbers. Note that the tem- 
perature contour distortion from the conduction profile increases with 
Rayleigh number as expected. The higher Rayleigh number temperature 
contours show a strong distortion from the low Rayleigh number results. The 
streamline contours, however, are basically unchanged over the range of 
Rayleigh numbers, except for a relatively slight asymmetry introduced at 
the higher Rayleigh numbers. The flow pattern consists of a single con- 
vection cell for all of the Rayleigh numbers shown. 
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The variations in temperature, streamline and velocity contours with 
time are shown in Figs. 6, 7 and 8 for Ra s 10,000. As expected, the dis- 
tortion in the temperature field increases with time toward the steady stale 
condition. The single convection cell is also shown to shift gradually from 
the symmetry characteristic of low Rayleigh numbers toward the asymmetry 
characteristic of high Rayleigh numbers. 
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JONCLUSIONS 

There are a number of mechanisms other than gravity- induced con- 
vection affecting the quality of crystal growth. The results obtained in 
this study, therefore, must be considered in conjunction with other infor- 
mation in evaluating the effect of gravity on the outcome of the Lal/Krocs 

experiment. The results of this study indicate gravity induced convection 

-5 -3 

velocities ranging from approximately 2x10 to 2x10 cm/sec for 
typical orbital gravity loads of 10"^ to 10"^ g e> Although these velocities 
appear to be extremely small, they are comparable to the expected crystal 
growth rate velocities and molecular diffusion velocities. 
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Part 2 

NUMERICAL SIMULATION OF NATURAL CON- 
VECTION WITH GRAVITY -SHIFT IN CIRCULAR 
CYLINDERS IN LOW GRAVITY 
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ABSTRACT 

A numerical analysis was performed to simulate natural convection in 
circular cylinder enclosures with gravity shifts similar to what might be ex- 
pected on board an orbiting spacecraft. The particular cases investigated 
were for a vertical Rayleigh number (based on the vertical component of 
gravity) of 1000 with superimposed horizontal Rayleigh numbers of +1000 to 
^ 000. The initial temperature gradient was taken to be in the direction of 
the horizontal. For a given set of vertical and horizontal gravity components, 
the results at steady state were found to be the same regardless of the order 
in which the two components were imposed. A supercritical horizontal 
Rayleigh number case with no vertical component was simulated to investi- 
gate transient development of convective flow from rest. The convective 
flow was found to grow exponentially from a small perturbation until steady 
state is approached. The numerical results for both subcritical and super- 
critical Rayleigh numbers were found to agree closely with results predicted 
by a simple analytical model derived for simultaneous loading of both vertical 
and horizontal components. 
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v 
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NOMENCLATURE 

Description 

cylinder diameter s 2 R 
gravity force 
horizontal g component 

vertical g component 


radial distance 
cylinder radius 

Rayleigh number 


r 3 ATd‘ 
va 


horizontal Rayleigh number 
vertical Rayleigh number 


temperature 

initial mid-point temperature 

temperature difference across circular 

cylinder 

time 

2 

dimensionless time = pt/R 
velocity 

spatial maximum velocity 

v max * or zero horizontal gravity component 


v - v 
max max 

o 

I 2g 

dimensionless velocity v * v 

g 3 AT a 

rectangular coordinates (Fig.l) 
thermal diffusivity 

volumetric coefficient of thermal expansion 
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Symbol 

0 

P 

V 

p 


NOMENCLATURE (Concluded) 

Description 

polar angle (Fig.l) 

dynamic viscosity 

kinematic viscosity, p/p 

stream function 

density 

p at T = T 
r o 
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INTRODUCTION 

Some of the fluid mechanics experiments planned to be performed 
aboard the forthcoming Space Shuttlc/Spaeelab flights are intended to investi- 
gate the use of the near -zero gravity conditions in orbit for various materials 
processing applications. The near elimination of gravity -induced convective 
stirring under orbital conditions is expected to yield various high quality 
specialized products that cannot easily be produced under terrestrial gravity 
conditions 

A number of theoretical and experimental studies have been made of 
convective flows within various geometric enclosures. Batchelor (Ref. 1) 
and Weinbaum (Ref. 2) investigated, respectively, rectangular and circular 
two-dimensional enclosures using non-numerical analytical methods for 
steady state conditions. Dressier (Ref. 3) recently extended Weinbaum' s 
steady state results for the cylinder to include a transient solution. These 
analytical results are valid only for low Rayleigh numbers. Robertson and 
Spradlcy (Ref. 4) performed numerical analyses to show that the low Rayleigh 
number analytical theory is valid for Rayleigh numbers up to approximately 
1000. All of these investigations have assumed terrestrial-like constant 
gravity conditions,, 

In the near -zero gravity environment aboard orbiting spacecraft, the 
small accelerations that simulate residual gravity conditions are likely to 
change with time in both magnitude and direction. The purpose of this numer- 
ical study is to predict convective flows in circular enclosures with gravita- 
tional loads that do not remain fixed in time. The Lockheed-developed General 
Intcrpolant Method (GIM) computer code (Ref. 5) was used in the numerical 
computations. The computations were performed on the NASA -Langley Cyber 
203 system. 
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PROBLEM FORMULATION AND NUMERICAL SIMULATION 

The problem investigated is that of two-dimensional natural convection 
within the circular cylinder enclosure shown in Fig. 1. 




Fig. 1 - Geometry for Circular Cylinder Enclosure 


The initial temperature distribution was based on a linear horizontal gradient 
in the positive x direction, with the boundary points held constant in time: 


T(r,o) = T q + (AT/2)x/R 
T(R,t) = T q + (AT/2) cose 


( 1 ) 


I 

I 

1 


The gravitation loading consisted of two components, a vertical component 
in the negative y direction and a horizontal component in either the positive 
or negative x direction. We investigated applying the gravitational components 
in the following sequences: (1) the vertical component applied until steady 
state, followed by the horizontal component superimposed to produce a new 

2 
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steady state; (2) both components applied simultaneously until steady state; 
and (3) the horizontal component applied first followed by the vertical com- 
ponent. 


The numerical simulation was based on a nodal point distribution 
generated by the GIM code geometry module by specifying an array of 20x20 
area elements with 21x21 nodal points. The circle was treated as a four- 
sided figure, each side being a quarter -circle arc. Using stretching functions, 
the circular area was divided into generalized quadrilateral elements with 
curvilinear sides, with the nodal points located at the four corners of each 
element. The circular geometry is shown in Fig. 2 with the computational 
grid network superimposed. 

For convenience in the numerical simulation, we assumed a cylinder 
radius R of 1 cm and a temperature difference AT of 100 C. The gravity 
components were changed accordingly to yield the proper Rayleigh number 
values. The fluid was assumed to have the thermophysical properties of 
water and to behave as a Boussinesq fluid in its thermal expansion character- 
istics. The thermophysical properties used in the numerical simulation are 
listed as follows: 


Property 

Value 

Viscosity, p. 

Thermal conductivity, k 
Density, p 
Specific heat, 

Thermal expansion coefficient, 0 

1 centipoise 
0.00143 cal/cm-sec-C 
1 gm/ cu cm 
1 cal/gm-C 

2.07 x 10' 3 4 /C 
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RESULTS 

Results were obtained Cor a vertical Rayleigh number, Ra^, o£ 1000 and 
and horizontal Rayleigh number s, Ra^t of +1000, +2000, and _+ 5000. We found 
that the final steady state result for a given set of vertical and horizontal 
gravity components is the same regardless of the order in which the two 
components are applied. The steady state spatial maximum velocities for 
the various Rayleigh number combinations are presented in Fig, 3 in terms 
of a percentage deviation from the Ra^ -0 case. These results are com- 
pared to a simplified analytical theory derived in the appendix for the case 
of simultaneous impulse applications of the vertical and horizontal gravity 
components. Reasonable agreement with the analytical prediction is shown 
in this comparison, particularly for low horizontal Rayleigh numbers. The 
numerical results show an increase in the spatial maximum of about 150% 
for a positive horizontal Rayleigh number of 5000 and a decrease of about 
40% for a negative Rayleigh number of 5000. Note that the percent deviation 
depends only on the magnitude of the horizontal Rayleigh number. 

The simplified analytical theory, in agreement with Weinbaum (Ref. 2), 
predicts that, for a zero vertical Rayleigh number, convective flow will not 
develop for horizontal Rayleigh numbers less than 9216. The GIM code, 
however, probably due to numerical noise in the program, will generate a 
small convective flow at subcritical Rayleigh numbers. Superimposing a 
vertical gravity component on these results leads to a duplication of the 
results obtained by applying the gravity components either in the reverse 
order or simultaneously. 

The temperature contours at steady state are presented In Fig, 4 for 
a vertical Rayleigh number of 1000 and horizontal Rayleigh numbers of 0, 
+1000 and 4-5000, The positive horizontal components arc shown to noticeably 
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increase the distortion in the isotherms as expected. The negative horizontal 
components decrease the distortion somewhat, but not as noticeably as in the 
positive cases, 

The velocity contours at steady state are shown in Fig. 5, again for a 
vertical Rayleigh number of 1000 and horizontal Rayleigh numbers of 0, 

+ 1000 and +5000. The exhibited flow appears to be nearly perfectly circular 
in all cases, with the maximum velocity indicated roughly near the r « 0.5 to 
0,6 radial position. The theoretical maximum occurs at r = l/^3 s= 0,58 
(see Appendix). Some of the departure from circular symmetry in the velocity 
contours is probably due to the numerical noise in computing such small ve- 
locities. The circularity of the flow is more clearly indicated by the steam- 
line contours shown in Fig, 6. A nearly perfectly circular flow pattern is 
shown in all cases. 

A supercritical horizontal Rayleigh ni mber case was simulated to 
investigate the transient development of convective flow from rest. The 
Ra^j = 10,000 case was selected with no vertical component, A small ini- 
tial velocity perturbation in the flow field is required to start the develop- 
ment of convective flow, According to the simplified theory outlined in the 
Appendix, this velocity perturbation should grow exponentially until the sim- 
plified theory is no longer applicable. We started the numerical simulation 
with a velocity perturbation of 10“ 10 cm/sec in a counterclockwise direction. 
The resulting flow development is shown in the plot of spatial maximum ve- 
locity as a function of time in Fig. 7. In the first stage of the numerical sim- 
ulation, the computed velocities rapidly increased to levels beyond that 
expected from the simplified theory, then decreased to a lower level, At 
this point the flow began to rise exponentially approximately as predicted by 
the simplified theory. The first stage of the numerical simulation probably 
represents a numerical noise. After this period of adjustment, the simulated 
flow development takes on an apparently realistic trend. Note that the initial 
slope in the nearly linear portion of the logarithmic plot is 0.64 compared to 
0.89 predicted by the simplified theory for a Rayleigh number of 10,000. The 
flow peaks out and reaches a steady state dimensionless spatial maximum 
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Velocity Contours at Steady State for Various Horizontal Rayleigh Numbers 
)Ra„ = 1000, Velocity Units = 10~3 cm/s eel. 








2 


5 


6 


3 „ 4 

Time, t (dimensionless) 

Fig. 7 - Spatial Maximum Velocity as a Function of Time for Supercritical 
Horizontal Rayleigh Number 
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velocity of approximately 0,2 after a dimensionless time of approximately 6, 
This steady state velocity is of the same order as the 0,385 dimensionless 
velocity (Ref, 4) corresponding to a purely vertical gravity component. 
Theoretically, the time required to reach steady state depends on the magni- 
tude of the initial velocity perturbation, which is arbitrary. Our computed 
time to reach steady state, therefore, is not necessarily representative of 
actual supercritical cases involving the development of convective flow from 
rest, 


The development of convective flow from rest for the supercritical 
Ra^j = 10,000 case is illustrated by the contour plots of temperature and 
streamlines in Fig. 8. Steady state was reached in this case without breakup 
into multiple convective cells , 
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CONCLUSIONS 

For subcritical horizontal Rayleigh numbers, the imposition of a given 
set of vertical and horizontal gravity components yields the same steady state 
convective flow field regardless of the order in which the two components are 
imposed, For supercritical horizontal Rayleigh numbers without a vertical 
component, the convective velocity initially grows exponentially from a small 
perturbation until steady state conditions are approached, For both sub- 
critical and supercritical Rayleigh numbers, the computer simulation results 
agree reasonably well with predictions based on a simple analytical model. 
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Appendix 

All previous analyses, including Wcinbaum (Ref. 2), Dressier (Ref. 3 ) 
and Robertson and Spradley (Ref. 4) have demonstrated that the steady state 
flow pattern, is circular with the velocity given by 


Y~ v max t 1 " (r/W 2 ] (r/R) 


(A, 1) 


For simplicity, we assume that the transient flow maintains the same form as 
Eq. (A.l) with v varying with time. 

Based on the shape of the temperature isotherms noted in Ref. 2, we 
assume that the temperature distribution is given by 


T s T q + ^sjr- (r/R) cosQ T' 

where T' is a perturbation term due to convective flow given by 


(A. 2) 


X'sA [ 1 - (r/R) 2 ] (r/R) sinO 
“A [1 - (x 2 + y 2 )/R 2 ](y/R) 


(A. 3) 


Along the y axis (0 ” tt/ 2), the perturbed temperature distribution defined by 
Eq. (A. 3) follows the same form as the velocity distribution given by Eq. (A.l). 
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Wo determined the constant, A, in Eq. (A. 3) by satisfying the steady state 
conduction equation shown below. 


v . VT a a V 2 T 


(A. 4) 


at the point of maximum velocity and temperature perturbation along the y 

|a»as*» ’ 

axis , r - 1/y 3 (x = 0, y = l/\3 ). A more accurate method of accomplishing 
this perhaps, would be by minimizing error in Eq. (A. 4) over the entire circular 
region, rather than concentrating on the maximum velocity and temperature 
perturbation point. For the simplified solution that we arc seeking, however, 
the approach we selected seems appropriate., Finding the temperature gradient 
and Laplacian at x = 0 and y =; 1/^2 fron't Eqs. (A. 2) and (A, 3) and solving 
Eq. (A. 4) yields: 


^3 v max R AT 


(A.5) 


The torque, My, on the fluid mass due to viscous shear stress on the 
interior cylinder surface is given by 


My - Zlt m 


r = R 


&7T a[T II V R 
^ f max 


(A. 6) 


We take the counter-clockwise direction as positive, since the convective flow 
will assume this direction. 

The fluid density distribution is based on a Boussinesq variation with 
temperature: 

p = P o [1 - P(T - T o )] (A.?) 
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which, from Eq. (A. 2), becomes: 


P ~ P Q jl - 3 [^jf (r/R) cosG + T'j J 


The torque due. to gravity is found by: 




P (g v x + s H y) «is 


(A. 8) 


(A *9) 


where the integration is carried out over the circular region. The vertical 
gravity component gy, is always positive and acts in the negative y direction, 
The horizontal component, gjj, may be either positive or negative according 
to its direction with respect to the x axis. Carrying out the integration of 
Eq. (A. 9) yields 



2 7T R 



P( g v c osG + g H sinO) 


r^ dr dO 


if 

8 


P 0 {3 AT R* 


(gy + 


V? 


v R 
max 


24 a 


g H ) 


(A. 10) 


The instantaneous angular momentum of the fluid mass is given by: 



r pv dS 


(A. 11 ) 
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where, again, the integration is carried out over the circular region. Carry- 
ing out this integration yields 


2.71 R 



Pvr dr d0 



ttR p v „ 
*0 max 


(A. 12) 


By balancing the rate of change of angular momentum to the net torque 
on the fluid mass, a differential equation is obtained for the variation of v 
with time: 


dP/dt = M r + M 


g 1 T v 


^7 r r3 Pq Vax= §PoP AT r3 («V + 


V? 


g H 


24a max 


^ - 6 7T hr 


v (A. 13) 
max v ' 


Re-arranging Eq. (A. 13) yields; 


v = — ^4 Ra v - 24 
max 16 V3 R 3 V 


/, . Ra n\ 

^2 y 92 16 / v max 


(A, 14) 


which yields the solution 


v max = (v max> t = 0 ex P 


-24 




(a/R) Ra y 

384 ^ (1 - £&) 


1 - exp 


-24 1 


Ra ii\ 

WK) r 2\ 


(A. 15) 
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whore (v _ a „) is an initial small perturbation in v^_ . Note that for 
max^ _ q ‘ max 

zero vertical gravity component (Ra^i: 0), any initial small perturbation in 

'max w *^ damped f° r Rajj<9216. The value 9216, therefore, is the critical 

Rayleigh number for initiating convective flow for cases where the gravitational 

vector is in the direction of the temperature gradient. This value agrees with 

Weinbaum' s (Ref. 2) stability analysis for the horizontal cylinder (Weinbaum 1 s 

value of 576 is based on a definition of Rayleigh number which differs by a 

factor of 16 from our Rayleigh number). 


For Ray> 0 and Ra^j< 9216, Eq. (A. 15) yields the steady state result; 


max 


max 


1 - 


92 TE 


(A. 16) 


where v 


max 


o 


is the maximum velocity for a vertical component only (R^jj - 0); 


max. 


384 V3 


a 

R 


Ra 


V 


(A. 17) 


The fractional deviation in v due to a side load is given by; 

max ° J 


v - v 
max max 


max. 


Ra 


1 - 


H 
Ra H 

~ 92 TG 


(A. 18) 


For horizontal Rayleigh numbers approaching or greater than the critical 
value of 9216, this simplified analysis is no longer applicable. 
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Part 3 

EFFECT OF CONTAINER SHAPE ON NATURAL 
CONVECTION IN ENCLOSURES 
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Numerical computations were made of natural convection in two- 
dimensional enclosures of various shape but of equal cross sectional area. 
Results were obtained for enclosures with circular, half-circular and square 
shapes. For the shapes investigated, the maximum convective velocity was 
found to be about the same for all shapes, This suggests that, as a general 
rule, complicated container shapes can be approximated by simpler shapes 
of equal area for purposes of predicting the intensity of natural convection 
in experiment configurations. Additional computations are planned for 
other enclosure shapes, such as a triangle and rectangle. The complete 
results will be documented in the future as a Lockheed technical report 
and as a paper in the open scientific literature. 
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